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ABSTRACT The structural, electronic, and magnetic properties of Co,O (n = 2~~10) clusters 
have been systematically investigated within the framework of the generalized gradient 
approximation density functional theory. The results indicate that the O atom occupies the 
surface-capped position on Co,O (n = 2~10) clusters. The stabilities of the host clusters are 
improved by adding one O atom. Maximum peaks of the second-order difference energy of the 
ground-state Co,O clusters are found at n = 3, 6 and 8, indicating higher stability than their 
neighboring clusters. Compared with corresponding pure Co, clusters, the O-doped cobalt clusters 
have larger gaps between the HOMO and LUMO energy levels, indicating their higher chemical 
stabilities. In addition, the doping of O atom exhibits different influence on the magnetism of the 
clusters. This is also further investigated by the local magnetic moment, deformation charge 
density and partial local density of states analysis. 
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1 INTRODUCTION 


Transition metal oxide clusters are widely used in high temperature chemistry, materials 


science, microelectronics and nanotechnology areas, and have attracted extensive theoretical and 
experimental attention in the past ten years. For example, Clemmer et al.!!! reported experimental 
measurement of Sc-O bond energies and ionization potentials of ScO2 using guided ion beam 
reactions between ScO* and NO2. Wang er al.2! have calculated the structural and magnetic 
properties of small transition metal oxide clusters TM,Om (TM = Sc, Ti, V, Cr and Mn; n = 1,2; m 
= 1~6) by using ab initio density functional theory approach, and revealed that the geometries of 
clusters are closely correlated with the ratios of TM to O and the number of valence electrons of 


the TM atoms. Yang et al.l“! studied the structural and electronic properties of La,O (n = 2~12) 


clusters using the density functional theory (DFT), and found the doping of O atom prefers to stay 
outside the cluster, and can greatly improve the stability of small La, clusters. In addition, the 
structures, electronic, and magnetic properties of transition metal (M = Cr, Mn, Fe, and Ni) oxide 
clusters and their anions and cations have been studied largely!*!?!, Gutsev G. L. et al.!!3 studied 
the polarizability per atom of the MO monomers (M = Fe, Co, and Ni), and found that the doping 
resulted in a steep decrease of polarizability as compared to the atomic values. However, it is 
obvious that the existing theoretical study on Co,O clusters is not sufficient, and that the detailed 
structural and electronic properties of the large size of Co,O (n25) clusters are still unclear. 
Therefore, it is desirable to extend the study of Co,O clusters to larger size so as to glean a 
comprehensive understanding of their exotic physical and chemical properties. We hope that our 
work can provide powerful guidelines for the corresponding experiments and promote the 
functional design of the clusters cohesive material. In the past decades, on the experimental side, 
the mass selected cobalt oxide cluster anions Co,Om (n = 5~~21, m = 0-2) are studied by using 
photodetachment photoelectron spectroscopy (PES), indicating that O atom causes a minor 
influence on the electronic structures!'4!, Jacobson and Freiser!!>! have performed Fourier 
transform mass spectrometry investigation on the reaction of cationic Co dimers and trimers with 
O2. This was followed by a kinetic study of the reaction of small Co,* clusters (n = 2-9) with Oo. 
On the theoretical side, Kosuke et al.!!®! reported the results of structures of cobalt oxide cluster 
cations ComOn,* by ion mobility spectrometry (IMS), combined with the theoretical calculation in 


order to verify the structures of the clusters, and concluded that (CoO)3-5* ions are 


monocyclic-ring structures and (CoO)¢7* has compact tower structures. Therefore, structural 
transition from the ring to compact structures occurs at (CoO) *. The ground states of CoO, (n = 
1~4) and their anions [CoO,]" were studied by density functional theory, with the BLYP 
exchange-correlation functionall!”]. However, the researches on single O atom doped cobalt 
clusters are seldom investigated except Liu er al."'®! who study the structure and magnetic 
properties of the small-sized Co,O (n = 1~5) clusters by using the density functional approach. 
As far as we know, there is no systematical first principle calculation on the large size of 


Co,O (n = 2~10) clusters. Our work aims to make clear how the O atom binds with the Co, (n = 
2~10) clusters by theoretical calculation with density functional theory method, and try to 


explore how the stability changes after O atom’s doping with cobalt clusters, and how the 
magnetic properties of the host clusters are influenced by doping, finally to disclose the physical 
origin of the magnetic behavior of Co,O clusters. We hope that our research work would provide a 
certain theoretical guidance for the practical application of cobalt oxide clusters assembling 


materials. 


2 COMPUTATIONAL DETAILS 


All calculations were performed at the DFT level with the DMol? package in the Materials 
Studio of Accelrys Inc!'?!. The exchange-correlation interaction was treated within the GGA using 
PW91 function (GGA-PW91)!°l. The double numerical basis set augmented with d-polarization 
and p-polarization functions (DNP) was utilized. For the numerical integration, a fine quality 
mesh size was used, and the real space cutoff of the atomic orbital was set at 5.5 A. The 
convergence criteria for structure optimization and energy calculations were set to fine with the 
tolerance for density convergence in SCF, energy, gradient and displacement 1 x 10º eV/A, 1.0 x 
105 a.u, 0.002 Hartree/A and 0.005 A, respectively. In geometry optimization procedure, we 
considered a number of initial structures including linear chains, planar and three-dimensional 
structures in this work to maximize our chance to find the ground state configurations of the alloy 


clusters. Meanwhile, we also look over other structures of La,O, ScnOm, YnO, Mn,O, Nb,O, 


Co,Rh, Co„Mn, ConV and Co,Fe®: 21-28] clusters. In addition, owing to the spin polarization, all 
optimized geometry was optimized again at various possible spin multiplicities. Consequently, the 
number of possible initial isomers increased very rapidly with the size increase of clusters. The 
calculations are implemented until the minimum energy is reached. We confirm the stability of the 
lowest-energy structures as minima of the potential energy surface by considering no imaginary 
frequency. 

To check the validity of the computational method in our work, we first perform the 
calculation on the Coz and CoO dimers using PW91, PBE and BLYP, respectively. As listed in 
Table 1, we choose the PW91 method finally because the results are consistent with previous 
theoretical and experimental data well, and these data are sufficient to show that it is effective and 


reliable to study Co,O (n = 2~10) clusters with this method. 


3 RESULTS AND DISCUSSION 


3.1 Geometrical structures of the Co,O (n = 2~10) clusters 

The low-lying geometries, symmetries, magnetic moment and energy relative (AE) of Co,O 
(n = 2~10) clusters are shown in Figs. 1 and 2. The lowest energy structures of pure Co, clusters 
are also included for comparison. The lowest frequencies of Co, and Co,O clusters are listed in 
Table 2. 

For Co,O clusters, a number of different stable configurations for each size cluster are found, 
but only three of them are shown in this paper. In the case of Co2O cluster, the lowest-energy 
structure is an isosceles triangle (C>,) with the O atom located at the apex. It has a total magnetic 
moment of 4 ug, which is consistent with previous theoretical calculation results exactly!!?!, A less 
stable structure with 0.106 eV higher energy than the most stable one has more relaxed Co—Co 
and Co-O bonds to bring with a higher magnetic moment of 6 us. 

The most stable Co30 cluster adopts a planar structure (C2) with 7 ug of total magnetic 
moment. Its metastable structure (3 us), which shares the same geometric structure with the 


lowest-energy of Co30, is energetically less favorable than the most stable one by just 0.022 eV, 


so they are two energetically degenerate states and magnetic bistable states. The third stable 
isomer is a tetrahedron (Cs) which is different from the ground state and metastable state totally. It 
is 5 ug in magnetic moment and 0.107 eV less stable in energy than the lowest energy structure. 

A triangular bipyramid (8 us) with the O atom located at the apex is found to be the most 
stable structure of Co4O with Cs symmetry, which is in good agreement with the previous 
theoretical results!!?!, In comparison, the second lowest-energy structure has similar geometry with 
the C3, symmetry and lies 0.021 eV higher in energy with respect to the ground state. The next 
low-lying energy isomer (C>,), which is energetically less stable than the most stable structure by 
0.034 eV, exhibits double coordinating of O atom to the tetrahedron. 

For the CosO cluster, we considered a lot of initial geometries including capped triangular 
bipyramid, octahedron and pentagonal pyramid. The lowest-energy structure is a capped triangular 
bipyramid (Cs) with an 11 ug of total magnetic moment, in which the O atom is located at the 
capped position. The second stable isomer (9 ug), which shares the same geometric structure and 
symmetry, is just 0.003 eV higher in energy than the lowest-energy structure. By contrast, the 
other stable isomer at a higher energy (0.012 eV) is different from the first two because the O 
atom tends to bind preferably to the bridge site of the CosO cluster. 

CosO cluster (14 ug) adopts an O-capped octahedron structure with C3, symmetry. In the 
metastable structure, the O atom coordinates with two connecting Co atoms and is coplanar with 
four Co atoms of the octahedron. The other low-lying structure (Cs) with 0.029 eV higher in 
energy shares the same configuration with the ground state, while has a smaller total magnetic 
moment of 12 ug which maybe results from the decrease in average bond length. 

The ground-state structure of Co7O (15 us) is a bi-capped octahedron with Cs symmetry, in 
which the O atom is located at the capped position. The next two low-lying states have the same 
structure (C2y) in which atom O lays coplanarly with the 5-membered ring of the hexagonal 
bipyramid. With total magnetic moments of 13 and 15 yup, their energies are higher than the 
ground-state structures by 0.044 and 0.051 eV, respectively. 

The lowest-energy structure of CogO cluster can be viewed as a tricapped octahedron (C;) 
with an 18 us of total magnetic moment and the O atom is capped on the Cos ground-state 


geometry. The metastable state is a bicapped pentagonal bipyramid (C>,) where O atom lies in the 


middle plane and is opposite to the bicapped position. Its energy is 0.021 eV higher than the 
lowest energy structure. A side-capped tetrahedral prism geometrical isomer (C;) with the O atom 
located at the capped position, which shares similar total magnetism with the ground state and the 
metastable state, is considered to be the next stable structure at a slightly higher energy (0.027 
eV). 

In terms of CooO, three low-lying structures are similar, which are all the bi-capped 
tetragonal antiprisms. The difference lies in the way how to cap. For the ground state structure (19 
uB) with Ca, symmetry, the bi-capped atoms locate oppositely. While the metastable state and the 
third state are the same completely on structure (C3,), in which the O impurity occupies an anterior 
capped position of pure Coo with notable instability (0.003 and 0.016 eV, respectively). 

In the case of Co100, the structure is different from the Cojo cluster, which tends to form large 
layered structure of Co clusters. In addition, three low-lying structures have the same 
configuration and symmetry. The metastable state (18 us) and the third state (14 up) are higher in 
energy than the ground-state structure (20 us) by 0.027 and 0.0055 eV, respectively. 

From the above analysis, geometry optimizations reveal that O atom tends to bind preferably 
to the bridge site of the Co»O and Co30 clusters. About geometries of Co,O (n = 5, 6, 7, and 8) 
clusters, we can see that O adds to bind preferably to three coordination sites with the ground state 
geometries of Co, (n = 5, 6, 7, and 8), forming surface-capped structures. Although O atom also 
binds with three Co atoms in clusters of Co40O and Co100, the way in which O atom binds is 
different. Their structures are based on the Co atoms rearrangement instead of the original 
ground-state structure of pure cobalt clusters. Additionally, for CooO, O binds to Co atom on the 
square edge in the most stable geometry. In a word, the results indicate that the O atom occupies 
the surface-capped position on Co, (n =2~10) clusters. Meanwhile, we also found that a number 
of similar low energy isomers exist as the low-lying structures of the Co,O (n = 2~10) clusters. 
3.2  Stabilities and electronic properties 

To investigate the influence of doping of O on the stability of cluster, we computed and 
compared the average binding energies per atom of the lowest energy structures of Co,O clusters 
and Co, clusters (Fig. 3). The atomic average binding energy E» (n) can be expressed by the 


following formulas: 


E»(Co,0) = [nE(Co) + E(O) — E(Co,0)|/(n + 1) (1) 


Ev(Con) = [nE(Co) — E(Con)\/n (2) 


where E is the total energy of the respective atom or cluster. n corresponds to the number of Co 
atoms in clusters. From Fig. 3, we can see that the average binding energies of Co,O clusters are 
larger than those of Co, clusters, indicating that the doping with O improves the stability of the 
main cluster. Moreover, the average binding energies of Co,O clusters are increasing with the size, 
which indicates that these clusters can continuously gain energy during their growth. Specifically, 
in the initial stage (n = 1~7), E» increases rapidly, while n goes from 7 to 10, the Ep exhibits 
obvious convergence trend. The change trend of Ep for pure Co clusters matches with the previous 
theoretical results well!?° 28, 291, 

In cluster physics, the second-order difference of cluster energies is another important quality 
that can reflect the relative stability of clusters. The second-order difference energy A2E(n) is also 


calculated using the following formulas below and plotted in Fig. 3: 


A2E(Con) = E(Con+1) + E(Con-1) — 2E(Con) (3) 


A2E(ConO) = E(Con+10) + E(Con10) — 2E(Con0) (4) 


where E is the total energy of cluster. As shown in Fig. 3, the A»E(n) exhibits an obvious even-odd 
oscillatory behavior for Co, clusters, indicating that the Co, clusters with even number of atoms 
are more stable than their neighboring sizes. A similar trend is observed in Co,O clusters except 
for Co30. The local peaks are found at n = 3, 6 and 8, which indicates that these clusters are 
relatively more stable than their neighboring clusters. Meanwhile, local minimum peaks appear at 
n = 5, 7 and 9 identically. Overall, the doping of O atom does not influence the change trend of 
the stability of clusters. 

In addition, the HOMO-LUMO gap is a characteristic quality of electronic structure and 
stability in clusters and is commonly used to measure the ability for clusters to undergo activated 
chemical reactions with small molecules. As shown in Table 2, the HOMO-LUMO gaps of Co,O 
are usually larger than those of Co, clusters except for n = 8, suggesting their high chemical 
inertness of the doped Co,O clusters. The reason for the larger Egap is that the O atom’s doping 


makes the Co-Co interaction weak in clusters, leading to the large Egap!*4!. That is consistent with 


the trend of the binding energies. Moreover, noticeable peaks of Egap are found at n = 3, 6 and 
10, indicating that these clusters must be more stable than the neighboring sizes. However, the 
stabilities of clusters for n = 5 and 8 are low. 

To understand the differences of HOMO-LUMO gaps between Co, and Co,O clusters, the 
pictorial representation of molecular orbitals for n = 2, 4, 6, and 8 is used and plotted in Fig. 4. By 


this way, the types of orbitals involved in the bond formation can be seen clearly. For Coz, the 
HOMO consists of an anti-bonding dz orbital along the bond axis, while the LUMO is also 
noticed to be an anti-bonding but from s orbital. With the doping of O atom, an obvious z bond 


formed by dz orbitals of Co atoms accounts for the HOMO of Co20 cluster, which results in 


energy reduction dramatically. In the LUMO of Co20 cluster, the dz? orbitals bond weakly with p 
orbital of O atom but still an anti-bonding orbital. Compared with Coo, the bonding in Co20 
reduces the energies of HOMO largely and LUMO slightly, so it is responsible for the largest gap 


difference among all the pairs (Co,/Co,0). The HOMO diagram of Coy is noticed to a strong bond 
between d», from para-Co atoms and s orbitals from the other two, exhibiting a distinct d-p 


hybridization characteristic. While, the LUMO locates 0.507 eV higher in energy and consists of a 
non-bonding from four d orbitals (d,. and di). When doped with one O atom, it exhibits 
nonbonding overlap in HOMO which is higher in energy than that of Cog. The LUMO of Co4O 
shows nonbonding as well as anti-bonding overlaps. As the cluster size n goes to 6, the HOMO 
consists of six dx being nonbonding and o bond is observed in LUMO, which contributes a small 
HOMO-LUMO gap. After doping with an O atom, a large z overlap is induced in HOMO and an 
in-plane z bond in LUMO. Meanwhile, the LUMO also has a non-bonding perpendicular to the z 
bond, resulting in a much lower HOMO and a higher LUMO. Both the HOMO and LUMO of Cog 
are mainly composed of nonbonding overlaps of d orbitals, while some obvious sideway overlaps 
between s, p and d orbitals appearing in HOMO, forming d-s, d-p and d-d delocalized o bonds. 
Moreover, there is an obvious x bond in the LUMO of CosO cluster. Therefore, the enhancement 
of bonding causes an overall reduction in energy level and a tiny decrease in gap compared to Cog 


cluster. 


We also perform Mulliken population analysis for the lowest-energy structures, and the 
atomic charges of O atom are presented in Table 2. It shows clearly that all of the O atoms are 
negatively charged. There is an obvious charge transfer from Co to O, indicating that the O atom 
is a Mulliken charge receiver. Furthermore, the charge transfer means the presence of ionic bonds 
between Co and O in Co,O clusters. It is an important reason for the higher stability of the doped 
Co,0O clusters! 

In order to understand the charge distribution between Co and O atoms of the doped clusters, 
we analyzed the charge difference density of Co30 and Co¢O clusters. The difference charge 
density is defined as the difference between the total charge density of the cluster and a 
superposition of atomic charge densities with the same spatial coordinates in the cluster. As is 
shown in Fig. 5, the red color means charge accumulation, while the blue represents charge 
depletion. Strong charge accumulation is observed around the O atom both in Co30 and Co¢O, 
while the depletion of charge is found from vacuum to Co site. Significantly, the electron 
accumulation is highly localized, signifying the characteristic of ionic bond between O and its 
neighboring Co atoms. Additionally, we can conclude that the extent of charge accumulation on O 
atom is greater in CogO than in Co30. 


3.3 Magnetic properties 


Following, we will focus on the magnetic properties of Co,O (n = 2~~10) clusters. The total 


magnetic moment of the ground-state Co,O (n = 210) clusters has been calculated and the 


results are presented in Figs. 1 and 2. For comparison, the values of pure cobalt clusters are also 
plotted. Clearly, the total magnetic moments of Co,O clusters increase with the size, which shows 
a similar trend with Co, clusters. For n = 2, 3, 6, 7 and 10, the total magnetic moments of clusters 
even remain constant after doping with the O impurity, and those of the Co, (n = 4, 5) clusters 
decrease while the values increase for n = 8 and 9. Interestingly, whether the moment increases or 
decreases, it keeps a difference of 2 up. The magnetic behavior of Co,O (n = 1, 3, 5) clusters is 
similar to that of the previous theoretical work"'®). 

In order to explore the magnetic origin further, we selected the ground state Co,O (n = 5, 


7, 9) clusters as representatives to discuss the influence of the doping of O atom on local 


magnetic moment of the main clusters (Fig. 6). For comparison, the local magnetic moment 
per atom of corresponding pure Co, clusters is shown in Fig. 6. As is evident from the figure, 
the magnetic moments of the Co,O clusters are chiefly derived from the Co atom and the 
contribution of O atom is very small. Additionally, the magnetic moment reduction or 
enhancement upon doping stems from the reduced or increased local magnetic moments on 
the Co atoms. For instance, as for the CosO cluster, the local magnetic moment of Co atom is 
reduced because of the doping of O atom, resulting in decrease on the total spin moment in 
contrast with the corresponding pure Cos cluster. In the case of n = 9, accompanying the 
enhanced local magnetic moments on Co in the doped cluster, the total magnetic moment 
increases. For the cases of n = 7, the O’s doping redistributes the spin moment on the Co 
atoms and the local moments decrease for the Co atoms near O and increase for those at a 
distance overall, thus the total magnetic moments of the clusters remain unaltered before and 
after loading O. A similar phenomenon also appears in previous study on the magnetism of 


atomic oxygen adsorbed Sc, (n = 2~14) clusters! 
To further understand the magnetic mechanism of Co,O (n = 1~~9) clusters, we provide the 


partial (PDOS) and local (LDOS) contribution of different orbital components of CosO, Co7O and 
CooO clusters as representatives. As is evident from Fig. 7, firstly, we note that the magnetism of 
Co,O is mainly from the contribution of Co atoms in which their electronic states below Fermi 
level come primarily from d electron state and the contributions from s and p electron states are 
very small. For the O atom, below the Fermi level, the integral area of spin-up p electrons is 
nearly equal to the spin-down p electrons in LDOS, indicating that the doping of O atom can 
hardly affect the magnetism of clusters, which coincided with the above analysis about local 
magnetic moment of clusters. In addition, between the Fermi level of —4 and —6 eV, there is a 
clear indication of hybridization between O 2p and Co 3d, resulting in the rearrange of charge 
density and in turn influences the magnetic moment of cobalt clusters. Generally speaking, the 
relative shift between the spin-up and spin-down bands can indicate the degree of spin exchange 
splitting. Moreover, the greater spin exchange splitting, the larger spin polarization and magnetic 
moment of cluster”, As shown in Fig. 7, for CosO cluster, we find that the spin-up d electrons 


integral area of CosO/Cos is evidently smaller than that of the pure Cos cluster. Meanwhile, the d 


electron peak of spin-down of CosO cluster is broadened around —5.5 eV below the Fermi level, 
which further counteracts the spin-up d electrons, indicating that the spin exchange splitting of 
DOS is decreased obviously. So, the total magnetic moment is decreased after O doping the pure 
Cos cluster. Compared with Co7 cluster, there is no much change on spin exchange splitting of 
Co70/Cos, and the differences of integral areas between spin-up and spin-down d electrons in pure 
Co; cluster and Co70/Co; are similar, so the total spin moment of Co7O cluster remains 
unchanged. In addition, we find that the spin exchange splitting of CooO/Coo is obvious, and the 
integral area of spin-up d electrons is evidently larger than that of spin-down below the Fermi 


level in LDOS, resulting in the increasing magnetism of the doping-O cobalt clusters. 


4 CONCLUSION 


By using the first-principles DFT-GGA calculations, the geometries, stabilities, and magnetic 
properties of the O doped Co, clusters have been systematically studied, with the results 
summarized as follows: 

(1) In the most stable structures of Co,O clusters, the O atom prefers to occupy the 


surface-capped position of Co, (n = 2~~10) clusters. Geometry optimizations reveal that O atom 


tends to bind the bridge site in the Co2O and Co30 clusters; the O atoms are capped by three 
coordination in Co,O (n = 5, 6, 7 and 8) clusters; for CooO cluster, the O atom is capped by four 
coordination, in which the corresponding Co, maintains its original ground state structure of pure 
cobalt cluster. In addition, by analyzing the most stable geometries of Co4O and Co100 clusters, 
we also found that the O atom prefers to bind with Co atoms in a way of three-coordinate, but 
Co atoms appear obvious rearrangement, and the Co100 cluster exhibits a layer structure 
model. 

(2) From the analysis of the second-order energy difference, we can conclude that Co,O 
clusters at n = 3, 6, and 8 possess relatively higher stabilities than their neighbors, which means 
that the doping of O atom improves the stability of the Co30, CosO and CogO clusters. 


(3) By analyzing the HOMO and LUMO of Co, and Co,O (n = 2, 4, 6, 8) clusters, we can 


conclude that the gap is mainly determined by the bonding strength between atomic orbitals. 

(4) Compared to those of pure cobalt clusters, the magnetism calculations show that the total 
magnetic moments for Co,O clusters emerge different changing trends. Because of the 
hybridization between the 3d state of Co and the 2p state of doping-O, the local magnetic 
moments of Co atoms changed. However, the contribution of spin moments from O atom is 
negligible in Co,O clusters. Thus, the magnetic moment reduction or enhancement by doping 


can be attributed to the decrease and increase of the local magnetic moments of the Co atoms. 
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Table 1. Calculated Bond Lengths and Binding Energy per Atom of the Coz 
and CoO Clusters Using Different Density Functional Theory Methods 


Co, CoO 
Methods E E 

R(A) Ey (eV/atom) R(A) Ep (eV/atom) 

PBE 2.139 2.791 1.642 3.453 

PW91 2.134 2.851 1.642 3.489 

BLYP 2.155 2.626 1.655 3.306 

Theoretical 2.123/2.130128, 29] 2.7428) 1.59507] 

Experimental 2.3189 1.72101 1.63164 3.940.145% 


-SS 


Coz Dan dup CoO Ca 4p a Cay Ong 0.106eV 


Cos Dan lOp CoO C, Sup a Cy Sup 0.021eV b Cay Spg 0.034eV 


Cos Dyn 13up CoO C, Ipa aG Mtg 0.003eV b C, Hyp 0.012eV 
2.425 


Cos Op [41g CoO Cy, l4up a Coy l4up 0.004eV b C, 12up 0.029eV 
Fig. 1. Ground-state geometries of the corresponding bare Co, clusters and the lowest-energy 


structures and low-lying isomers of Co„0O (n = 2~6) clusters. Distances are given in A 


2.425 2.435 2.493 2.409 
1.776 1.794 
Co; Cy 15 ug CoO C, 15 ug a Cy 13 pg 0.004cV b Ca 15 pp 0.051 eV 
Cog Dry 16 pp CoO C, 18 ug a Cy 18 pg 0.021 eV b C, 18 ug 0.027 eV 
Coy Dyn 17 pip CoO Ca 19 ua a Cy, 19 pp 0.003 eV b Ca, 17 pp 0.016 eV 


Cojo Da 20 Hp CoO Ca 20up a Cy 18 up 0.027 eV b Cay 14 pp 0.055 eV 


Fig.2. The same as those in Fig. 1 but for n = 7, 8, 9 and 10. Distances are given in A 


Table 2. Frequencies (cm), HOMO-LUMO Gaps (eV) and Total Spin Magnetic Moments (us) of 


Con/ConO Clusters. The Charges (e) on O Atoms of ConO Clusters Are also Presented 


Cluster size (n) Frequency HOMO-LUMO gap Charge 
2 327.9/224.6 0.263/0.886 —0.681 
3 257.7/43.9 0.632/0.946 —0.680 
4 71.6/119.3 0.507/0.756 —0.693 
5 86.5/42.5 0.353/0.590 —0.694 
6 140.4/121.2 0.233/0.658 —0.695 
7 69.8/76.6 0.353/0.642 —0.691 
8 79.1/81.5 0.423/0.370 —0.689 
9 85.4/42.7 0.198/0.389 —0.683 
10 70.6/30.9 0.257/0.431 —0.692 


E, (eV/atom) 
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Fig. 3. Average binding energies per atom and the second-order difference energies 


of the most stable structures of Co, and Co,O clusters as a function of cluster size 
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Fig. 4. 3D isosurface HOMO and LUMO diagrams of Co, and CoO clusters (n = 2, 4, 6, 8) 


Fig.5. Calculated deformed charge density plots for Co3O (left) and CosO (right) 
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Fig. 6. Local atomic magnetic moments on each atom of the lowest energy structures of ConO (n = 5, 7, 9) 
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Fig.7. Partial (PDOS) and local (LDOS) density of states of Con and ConO (n = 5, 7, 9) clusters. 
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The broadening width parameter is chosen as 0.1 eV 
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The structural, electronic, and magnetic properties of Co,O (n = 2~10) clusters have been 
systematically investigated by density functional theory. The results indicate that the O atom 
occupies the surface-capped position on Co,O (n = 2~10) clusters. The doping of O atom on Con 
cluster improves the stability of the host cluster and Co30, Co6O and CogO clusters exhibit higher 
stabilities than their neighboring clusters. The charge distribution between O and Co in Co„O 
clusters shows that O atoms are of great ionic character. The magnetic moment reduction or 
enhancement by doping can be attributed to the decrease and increase of the local magnetic 
moments of the Co atoms. In other words, the doping-O atom almost has no effect on the total 
spin moment. The doping-O atoms contribute little to the total moments of Co,O clusters, but 
O atoms induce different influences on the local magnetic moments of neighbor Co atoms, which 


finally change the total moments differently. 
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